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Abstract 


Simultaneous  observations  of  an  active  region  located  near  central  meridian  were  obtained  with  the 
Very  Large  Array,  the  Solar  Mciximum  Mission  X-Ray  Pclychromator,  and  the  Beijing  Observatory 
magnetograph  on  18  December  1987,  during  the  Coronal  Magnetic  Structures  Observing  Campaign 
(CoMStOC).  An  asymmetric  loop-like  structure  connects  the  strong  leading  sunspot  with  a  nearby 
region  of  opposite  polarity.  Both  6  and  20  cm  emission  lie  along  this  structure,  rather  than  over 
the  sunspot,  with  higher  frequency  emission  originating  closer  to  the  footpoint  inside  the  sunspot. 
The  20  cm  emission  is  due  to  a  superposition  of  2nd  and  3rd  harmonic  gyroemission.  where  the 
field  strength  is  160  G  -  300  G,  while  the  6  cm  emission  is  due  to  3rd  harmonic  gv'roemission  from 
a  repon  where  the  magnetic  field  strength  ranges  from  547  G  -  583  G.  A  high  \-alue  of  the  Alfven 
speed,  ~  40,000  km /sec,  is  obtained  at  the  location  of  the  6  cm  source,  with  somewhat  lower  \’alues. 

10,000  -  20.000  km/sec,  at  the  location  of  the  20  cm  emission.  At  the  location  of  the  6  cm  source, 
the  plasma  temperature  diminishes  with  height  from  2.5  x  10®  K  at  ~  5000  km  to  1.3  x  10®  K  at 
~  15,000  km. 

The  X-ray  data  associated  with  an  area  of  trailing  plage  were  used  to  predict  the  brightness 
temperature  structure  due  to  thermal  bremsstrahlung  omission  in  the  6  and  20  cm  wavebands.  The 
predicted  6  cm  brightness  temperature  in  and  around  the  location  of  the  X-ray  peak  is  low.  consis¬ 
tent  with  our  lack  of  observed  6  cm  plage  emission.  The  predicted  20  cm  brightness  temperature  is 
consistent  with  that  observed  in  the  central  portions  of  the  plage,  but  the  high  20  cm  polarization 
requires  the  presence  of  cool  (Te  <  5  x  10’  K),  absorbing  plasma  overlying  the  hot  plasma  observed 
in  X-rays.  Two  different  models  for  interpreting  these  observations  arc  described. 
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1  Introduction 


Observations  of  coronal  loops  taken  with  the  Very  Large  Array  (VLA)  show  that  emission  at  dif¬ 
ferent  wavelengths  originates  at  different  heights,  with  longer  wavelengths  corresponding  to  higher 
levels.  Synthesis  maps  at  20  cm  wavelength,  for  example,  outline  loop-like  structures  that  stretch 
across  regions  of  opposite  magnetic  polarity  (Lang,  Willson,  and  Rayrole  1982)  while  observations 
at  shorter  wavelengths  (2  and  6  cm)  generally  indicate  sources  that  lie  along  the  legs  of  loops  (Alis- 
sandrakis  and  Kundu  1982;  Lang  and  Willson  1982;  Lang,  Willson  and  Gaizauskas  1983;  Strong, 
Alissandrakis  and  Kundu  1984).  The  20  cm  loops  are  similar  to  their  X-ray  counterparts,  for  the 
loops  detected  have  similar  sizes,  shapes,  and  brightness  temperatures. 

Although  the  X-ray  emission  is  attributed  to  thermal  bremsstrahlung  and  atomic  line  emission 
from  hot  (Te  ~  3  x  10®  K)  coronal  plasma,  the  microwave  emission  may  be  due  to  either  ther¬ 
mal  bremsstrahlung  or  to  gyroresonance  emi^ion  from  thermal  electrons  spiralling  along  coronal 
magnetic  field  lines,  in  order  to  determine  which  of  these  mechanisms  is  dominant,  it  is  neces¬ 
sary  to  know  the  electron  temperature  and  emission  measure  of  the  plasma  so  that  the  brightness 
temperature  at  a  particular  wavelength  can  be  calculated  and  compared  with  the  observed  values. 
Observations  of  soft  X-ray  spectral  lines  obtained  by  the  X-ray  Polychromator  (XRP)  aboard  the 
Solar  Maximum  Mission  (SMM)  satellite,  for  example,  can  yield  estimates  of  the  electron  tempera¬ 
ture.  emission  measure,  and  density  in  quiescent  loops.  These  parameters  have  been  combined  with 
plausible  estimates  for  the  magnetic  scale  height  and  the  harmonic  to  show  that  the  microwave- 
emitting  plasma  can  be  optically  thick  to  either  thermal  bremsstrahlung  or  gyroresonance  emission 
at  20  cm  wavelength  (Lang,  Willson,  and  Rayrole  1982;  McConnell  and  Kundu  1983:  Shevgaonkw 
and  Kundu  1984;  Gary  and  Hurford  1987). 

Direct  comparison  of  6  cm  and  soft  X-ray  emission  has  confirmed  the  importance  of  gyrorcso- 
nance  emission  at  this  wavelength.  However,  the  brightest  6  cm  sources  are  usually  not  associated 
with  the  brightest  X-ray  sources  and  the  dctadled  correspondence  between  the  radiation  at  these 
two  wavelengths  is  poor  (Schmahl  et  al.  1982:  Webb  et  al.  1983).  The  reason  for  this  is  that 
6  cm  gyroresonance  emission  typically  occurs  in  the  strong  magnetic  fields  above  sunspots,  wlule 
the  X-ray  emission  is  usually  most  intense  between  regions  of  opposite  magnetic  polarity.  Thermal 
bremsstrahlung  at  6  cm  is  often  too  optically  thin  to  be  delected. 

More  recent  comparisons  between  simultaneous  20  cm  and  XRP  images  have  shown  better 
agreement  between  the  emission  at  these  wavelengths,  although  there  are  significant  dilferenccs 
(Webb'efal.  1987;Xahg  et  al.  1987a,b:  Nitta  et  al.  1991  -  CoMStOC  I).  The  20  cm  maps  suggest 
-the, presence  of  both  cool  and  hot  plasmas  in  different  regions  that  are  not  detected  in  X-rays,  and 
there  is  X-ray  emitting  plasma  that  remains  invisible  at  20  cm.  Although  thermal  bremsstrahlung 
could  explaun  the  bulk  of  the  20  cm  emission,  the  detection  of  a  line-like  enhancement  in  the 
brightness  spectrum  at  the  apex  of  one  coronal  loop  suggests  the  dominance  of  gjTorcsonance 
emission  in  that  part  of  the  loop  (Willson  198-5). 

In  this  paper  we  describe  the  results  of  simultaneous  VL.-\  and  XRP  observations  of  a  solar 
active  region  obtained  during  the  Coronal  .Magnetic  Structures  Observing  Campaign  (Co.MStOC). 
The  primary  objective  of  CoMStOC  was  to  obtain  quantitative  information  about  the  plasma  and 
magnetic  field  structure  above  sdfar  active  regions  by  comparing  the  microwave  and  soft  X-ray 
data  with  theoretical  models  (Holman  and  Kundu  1985.  Brosius  and  Holman  1988.  1989)  and 
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phd^piienc  magnetic  field  measnrements.  In  Section  2  we  present  and  disdiss.  the.  coordinate 
obse^tions:  VLA  observations  at  four  different  firequendes  in  each  of  the  6  and  20  m  wavebands, 
XRP  spectral  and  ims^e  data,  and  a  photospheric  lon^tudinal  magnetogram  from  Bdpng  Obser¬ 
vatory.  In  Section  3  we  analyze  and  interpret  the  VLA,  XRP,  and  ms^etograph  data  in  terms  of 
physical  models.  In  Section  4  we  summarize  our  condusions  and  their  implications. 


2  Observations  and  Data  Analysis 

2.1  VLA  Observations 

ft. 

The  VLA  (B  configuration)  was  used  to  observe  the  active  re^on  AR  4906  (S34.5,  E17.7  at  1300 
UT)  between  1440  UT  and  2300  UT  on  18  December  1987.  The  entire  array  of  27  antennas  was 
used  at  four  frequendes  in  the  20  cm  band  (1373, 1420, 1550,  and  1665  MHz)  and  four  in  the  6  cm 
band  (4600, 4700, 4800,  and  4900  MHz)  with  bandwidths  of  12.5  MHz.  The  half-power  beamwidth 
of  the  individual  antennas  ranged  between  29^  and  34'.7  at  20  cm,  with  synthesized  beamwidths  of 
4".  At  6  cm,  the  primary  beamwidths  ranged  from  9'.7  to  10'.4,  with  synthesized  beamwidths  of 
P'.  One^pmr  of  frequendes  was  observed  for  about  5  minutes,  so  that  all  dght  could  be  observed 
in  about,  twenty  minutes.  The  data  were  sampled  every  10  s  and  calibrated  using  observations 
of  PKS1655+078.  Corrections  were  made  for  the  difference  in  '.he  signals  from  high  temperature 
noise ’diodes  located  on  four  of  the  antennas.  The  temperatures  of  these  sources  were  measured 
acrossieach  band  prior  to  the  observations.  The  data  were  edited  and  calibrated  with  the  standard 
solar  procedures  at  the  VLA,  and  used  to  make  synthesis  maps  of  both  the  total  intensity,  I,  and. 
Stokes  parameter,  V,  during  the  8.5  hour  interval  at  all  eight  frequencies.  Examination  of  the  fringe 
visilnDty  on  several  baselines  indicated  no  variations  that  could  be  attributed  to  solar  bursts  or 
interference. 

Figure  la  shows  the  frequency-averaged  20  cm  map  with  two  main  components:  the  source  on 
the  rightis  associated  with  the  sunspot  and  the  one  on  the  left,  which  is  ~  50%  left-hand  circularly 
pol.  rized,  is  associated  with  the  plage.  Individual  maps  at  the  four  20  cm  frequencies  (Fig.  2)  show 
thatthoe  is  a  systematic  difference  in  the  structure  of  the  20  cm  emission  at  these  frequencies.  The 
sharply  defined  peak  of  the  sunspot  source  systematically  shifts  by  ~  10"  in  the  direction  of  the 
6  cm  source  (Fig.  lb)  from  1375  MHz  to  1665  MHz.  This  is  a  manifestation  of  gyroresonance  as 
the  dominant  emission  mechanism:  as  the  observing  frequency  changes,  the  magnetic  fidd  strength 
required  to  stimulate  gyroresonance  emission  also  changes:  the  appropriate  fidd  strength  is  at  a 
slightly  different  location,  resulting  in  a  shift  of  the  radio  source  with  observing  frequency. 

In  the  B-configuration  VLA,  the  maximum  fringe  spadng  is  ~  225"  at  20  cm  and  66"  at  6 
cm.:  No  information  is  available  concerning  sources  larger  than  these  critical  limits,  so  possible 
extended,  low  brightness  temperature  structures  will  be  resolved  out  (see  CoMStOC  I  for  a  more 
detailed  discussion ).  The  6  and  20  cm  sources  seen  in  Figure  la.b  are  smaller  than  their  respective 
critical  sizes,  suggesting  that  no  emission  from  the  sources  seen  in  these  maps  was  resolved  out. 
Furthermore,  no  significant  changes  in  structure  or  brightness  temperature  (<  15%)  were  found 
between  the  maps  produced  using  the  standard  CLE.4N  algorithms  and  those  produced  using  the 
Maximum  Entropy  Method  (MEM)  (Cornwell  1986.  Bastian  1987).  The  MEM  algorithm  is  better 
able  to  restore  large-scale  structure  (as  long  as  it  is  smaller  than  the  maximum  fringe  spacing) 


than  is  CLEAN.  Thus  the  similarity  of  the  MEM  and  the  CLEAN  maps  indicates  that  there  is  no 
evidence  for  unresolved  large-scale  structure  in  either  the  6  or  20  cm  maps.  Neither  CLEAN  nor 
MEM,  of  course,  is  capable  of  restoring  the  quiet  sun  background  contribution,  which  is  -n*  4  x  10^ 
K  at  6  cm  and  ~  1  x  10®  K  at  20  cm. 


2.2  XRP  Observations 

The  XBP  Flat  Crystal  Spectrometer  (Acton  et  al.  1980),  collimated  to  14"  FWHM,  could  build 
images  by  rastering  over  a  portion  of  a  4'  x  4'  area  with  a  minimum  pixel  spacing  of  5",  or  scan 
a  series  of  spectral  lines  at  one  or  more  spatial  locations.  XRP  mapped  AR  4906  in  the  Fe  X\TI 
resonance  line  at  15.01  A,  the  most  sensitive  line  at  active  region  temperatures.  This  map  (Fig. 
Ic)  shows  two  peaks  of  emission:  the  weaker  source  on  the  right  lies  between  the  sunspot  and  a 
nearby  region  of  opposite  polarity;  the  stronger  source  on  the  left  is  associated  with  the  trmling 
area  of  plage  (see  Fig.  Id).  ‘  ~ 

A  spectroscopic  scan  (between  13  and  19  .4)  was  done  at  the  brightest  Fe  X\'TI  pixel  (marked 
with  an  “X“  in  Fig.  1)  in  the  SMM  orbit  beginning  at  0343  UT  on  1987  December  18.  Analysis 
of  these  spectral  data  follows  the  procedure  outlined  in  CoMStOC  I.  The  emissi^ty  calculations 
from  Mewe,  Gronenschild,  and  van  den  Oord  (1985)  and  the  coronal  abundances  of  Meyer  (1985) 
were  used  to  calculate  the  electron  temperatures  from  five  selected  line  flux  ratios.  Elach  of  these 
ratios  varies  by  a  factor  of  10  or  more  in  the  temperature  range  of  interest  (2-4  xlO®  K);  this 
steep  gradient  makes  the  ratio  sensitive  to  small  changes  in  temperature  and  relatively  insensitive 
to  systematic  uncertainties  such  as  relative  flux  calibrations,  elemental  abundance  variations,  or 
errors  in  the  atomic  parameters.  (See  CoMStOC  I  for  a  detailed  discussion  of  these  calculations  and 
possible  sources  of  error.)  These  values  and  their  uncertainties  (statistical  only)  are  listed  in  Table 
1.  The  uncertainties  in  the  weighted  mean  include  systematic  as  well  as  statistical  uncertainties. 
The  column  emission  measure  corresponding  to  each  temperature  in  Table  1  was  computed  using 
the  0  Vni  line  flux,  which  has  an  emissivity  curve  with  a  broad,  flat  response  in  the  range  of  inter¬ 
est:  therefore,  llie  emission  measure  determined  using  this  line  is  relativeiy  insensitive  to  possible 
errors  in  temperature.  Once  again,  the  uncertainties  on  the  individual  entries  are  statistical  only 
(dominated  by  the  uncertainty  in  the  0  VUI  line  flux)  while  the  weighted  mean  includes  systematic 
as  well  as  statistical  uncertainties. 

XRP  mapped  AR  4906  at  the  peaks  of  five  bright  resonance  lines — 0  VIII,  Ne  EX,  Fe  XVII, 
Mg  XI.  and  Fe  XVIII — in  the  SMM  orbit  beginning  at  1134  UT  on  1987  December  18.  Ratios  of 
these  images  can  be  used  to  determine  the  spatial  structure  of  the  plasma  parameters  across  the 
region  (as  was  done  in  CoMStOC  I).  Unfortunately,  AR  4906  was  relatively  weak  and,  although 
Fe  XVII  was  observed  throughout  the  region,  significant  counts  were  observed  in  only  one  other 
line,  0  VIII.  near  the  brightest  central  areas  of  the  plage.  The  O  VIII/Fe  XVII  flux  ratio  varies  by 
only  a  factor  of  about  five  in  the  temperature  range  2-4  xlO®  K.  making  it  more  sensitive  to  errors 
in  flux  calibration,  elementtil  abundance,  and  atomic  parameters  chan  the  line  ratios  a\'aiiable 
at  the  brightest  Fe  XVII  pi.xel.  Even  though  the  data  were  taken  8  hours  apart,  the  resulting 
electron  temperatures  and  emission  measures  obtained  with  the  O  VIII/Fe  XVII  flux  ratio  from 
the  spectroscopic  and  imj^e  data  were  equal  at  the  same  location.  This  indicates  that  the  active 
region  was  stable  during  the  observation  period.  .Also.  GOES  soft  X-ray  data  show  no  significant 
activity  at  or  between  the  times  of  these  observations. 
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A  comparison  of  electron  temperatures  obtained  from  various  line  flux  ratios  at  the  location  of 
the  spectroscopic  scan  leads  to  the  conclusion  that  the  Fe  XVII  emissivity  is  inconsistent  with  that 
of  the  other  lines  by  a  factor  of  ~  2.75.  For  greater  detail  see  CoMStOC  I,  where  the  factor  was 
found  *  3  be  lower  (~  2.25),  but  the  difference  is  within  the  uncertainty  of  the  measurements.  This 
emissivity  “correction  factor”  is  used  to  obtain  temperature  and  emission  measure  maps  (Figures 
le  cind  If)  from  the  0  VIII/Fe  XVII  map  flux  ratios.  Maps  of  the  predicted  microwave  brightness 
temperatures  at  20  cm  and  6  cm  are  also  shown  (Figures  Ig  and  11.)  and  discussed  further  in 
Section  3.3.  These  maps  were  smoothed  with  a  3  x  3  boxcar  to  remove  statistical  noise,  and  the 
calculations  were  done  only  for  those  pixels  with  significant  count  rates,  >  3£7’,  in  both  lines.  It  is 
difficult  to  obtain  a  meaningful,  quantitative  estimate  of  the  error  in  the  map  temperatures  and 
emission  measures  other  than  to  note  that  they  are  expected  to  be  larger  than  at  the  Fe  XVH  peak 
(Table  1). 


2.3  Beijing  Magnetogram 

Figure  Id  shows  a  photospheric  longitudinal  magnetogram  of  .\R  4906  on  18  December  1987. 
obtrined  from  the  Beijing  Observatory  (courtesy  of  H.  Zhang  and  H.  Wang).  The  sunspot  is  the 
oval-shaped  pair  of  solid  contours  to  the  right  of  map  center.  The  maximum  line-of-sight  umbral 
magnetic  field  strength  is  ~2600  G,  and  the  noise  level  is  ~100  G.  Comparison  of  the  umbral  field 
strength  with  data  taken  at  the  Mt.  Wilson  Observatory  shows  good  agreement,  indicating  that 
the  magnetogram  did  not  saturate.  Just  left  (east)  of  the  sunspot  is  a  nearby  region  of  opposite 
polarity  with  a  maximum,  inward  directed  field  of  ~  790  G.  Farther  to  the  left  is  the  trailing  plage 
with  a  maximum,  inward  field  of  "-930  G. 

Figure  Id  also  shows  field  lines  obtained  with  the  potential  field  extrapolation  code  developed 
by  Sakurai  (1982);  this  code  solves  Laplace's  equation  to  calculate  current-free  magnetic  field  lines 
using  the  line-of-sight  photospheric  fields  as  boundary  conditions.  A  cluster  of  field  lines  in  the 
lower  left  portion  of  the  sunspot  terminates  in  and  around  the  nearby,  relatively  strong  region  of 
opposite  polarity.  Tiic.se  field  lines  outline  a  loop-like  structure  (or  bundle  of  loops)  which  is  quite 
a.syninietric.  For  several  representative  field  lines,  we  list  pairs  of  total  extrapolated  photospheric 
magnetic  field  strengths  (in  Gauss),  where  the  first  element  corresponds  to  the  footpoint  inside  the 
sunspot  and  the  second  to  the  footpoint  in  the  region  of  opposite  polarity:  (1800,33),  (1500,280), 
( 1400.550).  The  corresponding  photospheric  footpoint  separations  for  the  above  field  lines  art  5.06. 
4.64.  and  3.80  xlO"*  km. 


3  Results 

3.1  Analysis  of  the  Sunspot 

The  Sakurai  code  was  modified  to  obtain  the  total  extrapolated  magnetic  field  strength  in  a  grid 
of  locations  all  .at  the  .same  specified  height  above  the  photosphere.  Contour  maps  at  heights  of  0. 
5000.  10000,  and  15000  km  arc  shown  in  Figure  3.  The  vertical  gradient  of  the  total  magnetic  field 
strength  near  the  renter  of  the  umbra  for  heights  between  0  and  5000  km  is  ~  0.30  G/km.  (This 
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is  less  than  but  comparable  to  the  gradient  of  0.41  -  0.67  G/km  obtained  by  Henze  et  al.,  1982, 
for  the  longitudinal  mcignetic  field  in  the  transition  re^on  and  photosphere  of  a  sunspot  umbra.) 
The  average  horizontal  gradient  in  the  total  photospheiic  magnetic  field  strength  is  ~  0.11  G/km. 
Knowing  the  height  variation  of  the  maximum  total  extrapolated  sunspot  field  and  assuming  that 
this  field  can  be  described  by  that  of  a  vertically  oriented  point  dipole  buried  below  the  photosphere, 
we  obtain  a  depth  of  1.76  x  10"*  km  and  a  dipole  moment  of  7.45  x  10^®  G-km^. 

These  dipole  parameters  are  used  to  calculated  contours  of  constant  ms^etic  field  strength 
which  extend  well  into  the  corona.  Figure  4  shows  magnetic  field  contours  corresponding  to  the  1st 
through  5th  harmonics  of  1665  MHz  and  4900  MHz.  The  dipole  field  is  compared  with  the  field 
calculated  with  the  Sakurai  code.  Specifically,  the  radii  in  the  x-direction  of  several  field  strength 
contours  in  the  point  dipole  field  (Fig.  4)  were  compared  with  the  radii  of  these  same  contours  in 
the  four  extrapolated  maps  of  Figure  3.  At  these  four  heights,  the  corresponding  contour  radii  are 
equal  to  better  than  10%  for  magnetic  fields  ~  300  G  and  we  find  that  the  point  dipole  field  is  a 
good  approximation  of  the  extrapolated  potential  field  within  this  three-dimensional  structure;  this 
approximation  simplifies  the  calculation  of  microwave  emission  from  sunspot  models.  For  locations 
in  which  B<300  G,  nearby  pockets  of  ms^etic  field  perturb  the  sunspot  field  from  that  of  a  point 
dipole. 

Figure  5  shows  the  theoretical  I  and  V  maps  for  1665  MHz  amd  4900  MHz  obtained  using  the 
point  dipole  approximation  and  the  model  sunspot  atmosphere  described  in  Brosius  and  Holman 
(1989).  The  failure  of  this  simple  model  to  reproduce  the  sunspot  microwave  structure  provides 
insight  into  the  temperature,  density,  and  field  structure  above  the  sunspot:  1)  The  coronal  temper¬ 
ature  and/or  density  must  be  relatively  low  over  large  portions  of  the  sunspot  and  2)  the  microwave 
emission  must  emanate  from  some  coronal  structure  other  than  the  sunspot  in  order  to  explain  the 
observed  6  and  20  cm  I  and  V  structure.  The  loop-like  structure  is  explored  in  the  next  sub-section. 

An  overlay  of  the  intensity  maps  in  the  20  cm  waveband  (Fig.  2)  with  the  magnetogram  reveals 
that  the  20  cm  intensity  contours  do  not  cover  the  entire  sunspot.  The  potential  field  extrapolations 
(Fig.  3)  show  that  magnetic  fields  which  are  high  enough  to  produce  2nd  harmonic  gyroemission  at 
all  four  frequencies  in  the  20  cm  waveband  (245  -  297  G)  occur  at  heights  of  5000  -  15000  km,  where 
coronal  temperatures  (>  10®  K)  and  densities  (~  10®  cm”®)  are  expected  to  persist  (Papagiannis 
and  Kogut  1975,  Hildebrandt  et  al.  1987).  Indeed,  2nd  harmonic  gyroemission  is  expected  to 
emanate  from  the  portion  of  the  sunspot  for  which  the  observed  brightness  temperature  is  less  than 
the  lowest  observed  intensity  contour  (5.3  x  10®  K).  Ways  to  obtain  brightness  temperatures  this 
low  from  the  2nd  harmonic  gyroemission  mechanism  are  1)  the  density  in  the  emission  repon  is 
extremely  low  (<  10®  cm”®  for  a  2  x  10®  K  pla?ma),  2)  the  temperature  in  the  emission  repon  is 
less  than  5.3  x  10®  K,  or  3)  both  the  tempera ‘nre  and  the  density  in  the  emission  region  are  low. 
To  obtain  a  density  less. than  10®  cm”®  in  the  corona  over  a  large  fraction  (~  1/3)  of  a  sunspot 
area  seems  unlikely,  but  it  cannot  be  ruled  out.  Thus  we  find  that  a  lower  temperature  and/or 
density  exists  over  a  portion  of  both  the  sunspot  umbra  and  penumbra  on  the  side  of  the  sunspot 
away  from  the  center  of  AR  4906.  This  is  consistent  with  tite  results  of  Gary  and  Hurford  (1987) 
and  with  Strong,  .Alissandrakis,  and  Kundu  (1984). 


3.2  Analysis  of  the  Loop  Structure  Connecting  the  Sunspot  with  the  Nearby 
Region  of  Opposite  Polarity 

An  overlay  of  the  intensity  maps  in  the  20  cm  waveband  (Fig.  2)  with  the  magnetogram  field  lines 
(Fig.  Id)  reveals  that  the  peak  of  the  sunspot-assodated  source  in  all  four  of  the  20  cm  intensity 
maps  lies  along  the  loop  structure  joining  the  sunspot  with  the  nearby  region  of  opposite  polarity. 
Figure  6  shows  the  locations  of  the  centers  of  the  20  cm  peak  intensity  contours,  along  with  the  6 
cm  I  and  V  sources,  on  the  magnetogram.  The  higher  the  frequency,  the  farther  down  the  leg  of 
the  loop  (on  the  side  of  the  sunspot)  the  intensity  peak  occurs.  Assuming  the  emission  in  the  20 
cm  peaks  to  be  dominated  by  2nd  harmonic  gyroemission,  we  estimate  a  m^netic  fidd  gradient  ~ 
0.0068  G/km  along  the  loop  in  the  plane  of  the  sky.  Assuming  3rd  harmonic,  the  fidd  gradient  is  ~ 
0.0047  G/km.  Both  of  these  gradients  are  factors  of  2  to  3  less  than  the  gradients  calculated  along 
the  individual  ^trapolated  field  lines  which  pass  through  the  20  cm  emission  peaks.  Based  upon 
the  potential  fidd  extrapolation,  the  20  cm  peak  emission  re^on  is  comprised  of  a  superposition  of 
2nd  and  3rd  harmonic  gyroemission  arising  from  the  loop  structure.  The  magnetic  fidd  strength 
ranges  from  160  to  300  G  in  the  20  cm  source  re^on. 

The  6  cm  emission  originates  farther  down  the  leg  of  the  loop  (on  the  side  toward  the  sunspot) 
since  much  higher  fidds  are  required  for  gyroresonance  emission  at  6  cm  than  at  20  cm.  The 
polarization  structure  shows  two  V  sources  -  one  on  dther  side  of  the  I  source  -  such  that  the  I 
and  V  sources  ail  lie  roughly  along  the  length  of  the  loop.  The  polarization  is  '>'60%  dose  to  the 
edges  of  the  I  source.  Such  structure  appears  in  the  microwave  loop  models  of  Holman  and  Kundu 
(1985)  and  arises  from  the  gyroresonance  emission  mechanism:  the  ordinary  mode  (0-mode)  of 
a  given  harmonic  becomes  optically  thin  toward  the  harmonic  edge  more  r^idly  than  does  the 
extraordinary  mode  (X-mode).  Assuming  the  emission  in  the  6  cm  peadc  to  be  dominated  by 
3rd  harmonic  gyroemission  (see  next  paragraph),  we  estimate  am  average  magnetic  fidd  gradient 
(between  the  6  and  20  cm  source  regions)  >>'  0.020  G/km  along  the  loop  in  the  plane  of  the  sky. 
This  is  less  than  but  comparable  to  the  average  fidd  gradient  '>'  0.028  G/km  calculated  along  the 
individual  extrapolated  fidd  lines  which  pass  through  the  6  and  20  cm  emission  peaks. 

In  principle,  the  harmonic  number  of  the  6  cm  emission  can  be  determined  by  comparing  the 
observations  with  models.  Although  a  simple  dipole  loop  model  does  not  accuratdy  portray  the 
as3nmmetry  of  the  extrapolated  fidd  lines,  we  are  only  interested  in  reprodudng  the  emission  from 
one  leg  of  the  loop  (the  sunspot  side  of  the  loop  in  Fig.  6;  the  right  half  of  Figure  7),  so  the  magnetic 
fidd  structure  of  the  loop  elsewhere  is  of  little  consequence.  The  density  is  selected  so  that  the 
brightness  temperature  calculated  for  the  portion  of  the  loop  away  from  the  desired  6  cm  source 
remains  less  than  the  lowest  observed  6  cm  brightness  temperature  contour,  as  observed.  Both 
thermal  bremsstrahlung  and  thermal  gyroemission  were  included  in  the  calculations  and  theoretical 
I  and  V  maps  at  6  cm  were  produced  for  three  cases  where  the  2nd,  3rd,  and  4th  harmonic  emission 
region  was  in  the  footpoint  at  a  hdght  of  5000  km.  The  6  cm  I  and  V  structure  and  magnitude 
could  not  be  produced  with  4th  harmonic  emission.  For  2nd  harmonic,  the  density  required  is  low: 
'>'  5  X  10®  cm~^.  If  the  density  were  higher,  the  overlying  3rd  harmonic  emission  region  would  yield 
too  high  a  brightness  temperature  from  a  large  area  of  the  loop  and  would  destroy  the  desired  V 
structure.  For  3rd  harmonic,  the  density  is  '>'  10®  cm~®.  This  is  reasonable  in  light  of  the  sunspot 
atmosphere  model  of  Hildebramdt  et  al.  (1987)  where  a  density  of  •>'  8  x  10®  cm“®  occurs  5000 
km  above  the  photosphere.  The  emission  measure  for  the  model  loop  with  density  «>'  10®  cm~®  is 
'>'  10^"  cm~®,  less  than  but  comparable  to  the  value  of  ^  10^®  cm~®  estimated  from  the  very  weak 
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Fe  XVn  flux  measured  by  XRP  at  the  location  of  the  6  cm  source  (and  assuming  =  2.5  x  10® 
K).  The  corresponding  value  for  the  model  loop  with  density  ~  5  x  10®  cm~^  is  2  x  10^  cm~®, 
substantially  less  than  the  observed  value.  Thus  our  model,  combined  with  XRP  observations, 
indicates  that  the  observed  6  cm  I  and  V  structure  is  indeed  due  to  3rd  harmonic  gyroemission. 
Accounting  for  the  four  frequencies  at  which  the  6  cm  emission  was  observed,  this  translates  into 
magnetic  field  strengths  ran^ng  from  547  G  to  583  G.  At  the  location  of  the  6  cm  source  the 
potential  field  extrapolation  yidds  a  field  ~  600  G  at  a  height  5000  km,  consistent  with  our 
deduced  values. 

These  m^netic  field  and  density  results  can  be  used  to  estimate  the  Alfven  speed  and  its 
variation  within  the  loop.  The  magnetic  field  strength  required  for  the  highest  fireqnency  source 
in  the  6  cm  band,  583  G,  ^ves  an  dectron  gyro&equency  of  1.63  GHz.  The  6  cm  source  modd  is 
consistent  with  the  observations  when  the  dectron  density  is  between  5x10*  cm~^  and  2  x  10® 
cm“^,  giving  a  plasma  frequency  of  0.2  -  0.4  GHz.  Hence,  the  ratio  of  the  gyrofrequen<y  to  the 
plasma  frequency  is  found  to  be  in  the  range  4-8  and  the  Alfven  speed  falls  mthin  the  28,000 
-  57,000  km/sec  range.  (The  estimated  density  of  1  x  10®  cm“^  pves  an  Alfven  speed  of  40,000 
km/sec.)  The  160  G  -  300  G  range  for  the  magnetic  fidd  strength  in  the  20  cm  source  region  pves 
an  dectron  gyrofreqnency  ran^ng  from  0.45  GHz  to  0.84  GHz.  The  dectron  density  is  not  as  weQ 
established.  With  an  estimated  value  of  1  x  10®  cm~^,  the  ratio  of  the  gyrofreqnency  to  the  plasma 
frequency  ranges  from  1.5  to  3,  and  the  Alfven  speed  ranges  from  10,000  km/sec  to  20,000  km/sec. 
Hence,  the  .A.lfven  speed  in  this  coronal  loop  is  highest  near  the  sunspot  and  decreases  toward  the 
apex  of  the  loop. 

In  wave  heating  theories  of  the  corona,  the  heating  rate  is  proportional  to  the  Alfven  speed  and 
the  square  of  the  wave  vdodty  amplitude.  The  .Alfven  speeds  derived  here  are  significantly  greater 
than  the  generally  quoted  value  of  KKK)  •  2000  km/s.  This  implies  that  the  same  heating  rate  can 
be  achieved  with  a  smaUer  wave  amplitude.  However,  for  resonant  waves  in  a  loop  (as  in  resonant 
absorption  theoriw),  the  period  of  the  driving  wave  must  be  shorter.  For  a  loop  length  of  100,000 
km  and  an  -Alfven  speed  of  20,000  km/sec,  the  wave  period  must  be  10  sec  or  less.  Moreover, 
the  gj'rofrequency  is  generally  assumed  to  be  less  than  the  plasma  frequency,  the  opposite  of  what 
is  found  here  and  in  other  CoMStOC  results.  The  value  of  this  ratio  determines  which  plasma 
instabilities  are  important. 

-At  the  location  of  the  6  cm  peak,  where  the  brightness  temperature  is  ~  2.5  x  10®  K,  the  20  cm 
brightness  temperature  is  lower,  ~  1.3  x  10®  K.  Because  lower  frequencies  require  lower  fidds  for 
the  thermal  gj'roemission  mechanism  and,  because  the  magnetic  field  diminishes  with  increasing 
height  in  the  atmosphere,  the  20  cm  emission  most  orignate  higher  in  the  atmc^here  than  the  6 
cm  emission.  Indeed,  fields  high  enough  for  20  cm  2nd  harmonic  gyroemission  (245  -  297  G)  occur 
over  the  6  cm  source  at  heights  in  excess  of  15,000  km.  If  the  brightness  temperature  is  ~  1.3  x  1(^ 
K  and  the  2nd  harmonic  is  optically  thin,  the  20  cm  emission  should  be  polarized  since  the  0-mode 
is  less  optically  thick  than  the  X-mode.  This  polarization  is  not  observed  and  we  conclude  that  the 
observed  20  cm  brightness  temperature  is  due  to  optically  thick  2nd  harmonic  gyroemission  in  a 
region  where  Te  ~  1.3  x  10®  K.  This  means  that  the  plasma  temperature  at  the  height  of  formation 
of  the  6  cm  emission  is  higher  than  the  plasma  temperature  at  the  height  of  formation  of  the  20 
cm  emission.  Thus,  at  the  location  of  the  6  cm  intensity  peak,  the  plasma  temperature  diminishes 
ruith  height 
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3.3  Analysis  of  the  Plage 


The  electron  temperature  and  column  emission  measure  values  obtained  from  the  XRP  data  can  be 
used  to  compute  the  thermal  bremsstrahlung  microwave  emission  from  the  X-ray  emitting  plasma. 
The  optical  depth  due  to  thermal  bremsstrahlung  (free-free)  absorption  is 

A-.o  9.786  X  HAlT^ju) 

n^T}-^{u^VBCOs9Y  ’  ^  ^ 


where  T*  is  the  electron  temperature  in  K,  EM  is  the  emission  measure  in  cm~®,  v  is  the  obser\'ation 
frequency  in  GHz,  n„  is  the  index  of  refraction  ('v  l  for  plasma  and  field  conditions  here),  vb  is  the 
electron  gyrofrequency  {vb  —  2.8  x  where  B  is  in  Gauss  and  i/jg  is  in  GHz),  9  is  the  angle 

between  the  field  and  the  line  of  sight,  and  tne  minus  sign  refers  to  X-mode  while  the  plus  sign 
refers  to  0-mode.  The  thermal  bremsstrahlung  brightness  temperature  predicted  for  each  mode 


can  be  calculated  using 


(2) 


In  a  mj^etized  plasma,  the  optical  depth,  and  hence  the  predicted  brightness  temperature,  of  the 
X-mode  emission  is  expected  to  be  higher  than  that  of  the  0-mode.  The  predicted  total  intensity 
iPRED  jg  simply  the  average  brightness  temperature  of  the  two  modes  which,  depending  on  the 
direction  of  the  magnetic  field,  are  associated  with  right  circular  polarized  (RCP)  or  left  circular 
polarized  (LCP)  emission. 


At  the  location  of  the  Fe  XVII  peak,  where  the  XRP  spectroscopic  scan  data  used  to  generate 
Table  1  were  obtained,  the  electron  temperature  of  the  plasma  observed  by  XRP  was  2.95  ±  0.20 
xlO®  K  and  the  log  of  the  column  emission  measure  was  28.85  ±0.19.  These  values  entered  into  th^ 
equations  above  yield  the  predicted  intensities  listed  in  Table  2.  The  uncertainties  are  due  to  the 
errors  associated  with  the  electron  temperature  and  emission  mea,'ure  values  in  Table  1.  At  this 
same  location  in  the  plage,  the  VLA  observed  an  intensity  ~  1.2  x  10®  K  at  all  four  frequencies 
in  the  20  cm  waveband,  and  Ig  <  0.42  x  10®  K  in  the  6  cm  waveband.  Our  predicted  intensities 
in  the  G  cm  waveband  are  all  consistent  with  the  VL.V  observations.  Our  predicted  intensities  in 
the  20  cm  waveband  are  about  1.4  to  1.9  times  higher  than  the  corresponding  observed  intensities. 
Thus  it  appears  that  thermal  bremsstrahlung  plays  a  significant  role  in  the  observed  plage  emission. 

The  picture  becomes  considerably  more  complex  when  the  observed  polarization  is  considered. 
Table  2  lists  the  observed  intensity  and  the  LCP  and  RCP  brightness  temperatures  and 

Tg^^)  for  each  of  the  four  frequencies  in  the  20  cm  waveband  at  the  location  of  the  XRP  spectro¬ 
scopic  scan.  The  X-mode  corresponds  to  the  RCP  emission  for  outward  directed  fields,  and  to  the 
LCP  emission  for  inward  directed  fields;  the  0-mode  corresponds  to  the  LCP  emission  for  outward 
directed  fields  and  to  the  RCP  emission  for  inward  directed  fields.  The  LCP  emission  is  reasonably 
consistent  with  the  predicted  thermal  bremsstrahlung  brightness  temperature  (I^^^®).  However, 
the  observed  RCP  emission  is  considerably  lower,  by  a  factor  of  2  to  4,  than  the  predicted  value. 
Thus,  some  mechanism  is  required  to  explain  the  reduction  of  the  expected  RCP  emission  to  its 
observed  value.  The  only  way  to  reduce  the  expected  brightness  temperatures  of  ~  2  x  10®  K  to 
the  observed  value  ~  7  x  10®  K  is  to  have  cooler  absorbing  plasma  along  the  line  of  sight  between 
the  emission  region  and  the  observer. 


One  possible  model  which  explains  the  discrepancy  between  the  observed  and  the  predicted 
polarization  is  that  the  source  region  consists  of  two  layers  of  hot  plasma  separated  by  an  inter- 
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veiling  layer  of  cooler  plasma.  The  lowest  layer  contains  the  highest  column  emission  measure  and, 
therefore,  is  responsible  for  the  majority  of  both  the  X-ray  emission  and  the  predicted  thermed 
bremsstrahlung  microwave  emission.  The  intervening  layer  of  cooler  plasma  absorbs  out  the  mi¬ 
crowave  emission  from  this  lower  hot  layer,  leaving  only  the  microwave  emission  from  the  uppermost 
hot  layer  to  be  detected  by  the  VLA  (see  Figure  8a).  The  uppermost  hot  layer,  optically  thin  to 
thermal  bremsstrahlung,  is  embedded  in  a  magnetic  field  such  that  thermal  gyroemission,  capable 
of  producing  the  observed  high  polarization,  produces  the  observed  20  cm  microwave  emission.  The 
cool  intervening  layer  must  have  a  temperature  <  5  X  10®  K  in  order  to  reduce  the  predicted  RCP 
emission  at  1665  MHz  from  1.75  x  10®  K  to  its  observed  value  of  5  x  10®  K.  If,  for  example,  the  cool 
absorbing  plasma  has  a  temperature  of  10®  K,  it  needs  log  EM  =  27.5  to  get  r//  ~  5  at  1665  MHz. 
This  means  that  if  the  cool  absorbing  layer  forms  a  sheath  ~  3500  km  thick,  it  needs  a  reasonable 
density  ~  3  x  10®  cm~^  in  order  to  be  an  effective  absorber.  If  gyroabsorption  contributes  to  the 
opacity  of  this  cool  layer,  the  required  emission  measure  may  be  lower.  In  order  for  the  uppermost 
hot  layer  to  both  emit  strongly  in  the  20  cm  waveband  and  produce  polarization  as  high  as  50%, 
the  layer  must  contain  magnetic  fields  sufficiently  high  for  gyroemission.  From  our  potential  field 
extrapolation  at  and  near  the  location  of  the  XRP  spectroscopic  seem,  we  find  that  the  requisite 
magnetic  field  strengths  are  present,  indicating  that  the  observed  20  cm  emission  could  originate 
in  a  gyroresonance  emission  layer  between  heights  of  ~  5000  and  ~  15,000  km.  Since  the  radia¬ 
tion  is  emitted  primarily  in  the  X-mode  and  the  longitudinal  field  direction  is  inward,  the  emitted 
radiation  will  be  LCP,  as  observed  (Figure  2). 

Plasma  and  magnetic  field  parameters  in  the  uppermost  hot  layer  can  be  calculated.  The  optical 
depth  of  the  uppermost  hot  layer  can  be  estimated  by  solving 

Tg'^  =  Tc  exp{—Tgg)  -f  Te(l  -  exp(-Ty^)) 

for  Tgg,  where  T^^  are  the  observed  X-  (LCP)  and  0-  (RCP)  mode  brightness  temperatures 
at  a  particular  frequency,  Tc  is  the  brightness  temperature  of  the  optically  thick,  underlying  cool 
absorbing  plasma  (Tc  <  5  x  10®  K),  and  Tg  is  the  electron  temperature  of  the  uppermost  hot  layer 
(Te  >  2  X  10®  K).  Solving  for  we  obtain 

=  h[in-TMTc-T^''')]-  (4) 

Assuming  Te  =  3x  10®  K  and  Tc  =  1  x  10®  K,  we  obtain  =  0.80  and  =  0.15  for  the  observed 
brightness  temperatures  at  1665  MHz.  Now  the  optical  depth  due  to  gyroresonance  absorption  can 
be  written  (e.g.,  Kundu  1965,  Zheleznyakov  1970) 

=  0.052^^nei/-‘LB(1.77  x  10-'°re)’-Hl  ±  cos^)^  sin^*-^  (5) 

where  s  is  the  harmonic  number,  He  is  the  electron  density,  u  is  the  frequency,  Lb  is  the  magnetic 
scale  height  (B/yB),  Tg  is  the  electron  temperature,  0  is  the  angle  between  the  field  and  the  line 
of  sight,  and  the  plus  sign  refers  to  X-mode  while  the  minus  sign  refers  to  0-mode.  The  ratio 

tuh/'^uh  =  (1  +  cos0)V(l  -  cos9f  (6) 

yields  the  angle  between  the  field  and  the  line  of  sight  which  is  required  in  order  to  produce  the 
calculated  X-  and  0-mode  optical  depths  in  the  upper  hot  layer.  For  the  Ccise  discussed  above,  the 
calculated  angle  is  66°,  comparable  to  the  angle  obtained  from  the  extrapolated  coronal  magnetic 
field. 
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The  magnetic  scale  height  Lg  can  also  be  calculated  from  the  potential  field  extrapolation.  It 
is  0.94  X  10®  cm  in  the  2nd  harmonic  gyroemission  region,  1.4  x  10®  cm  in  the  3rd,  and  1.8  x  10® 
cm  in  the  4th.  Using  ~  0.8,  Te  =  3  X  10®  K,  i/  ~  1500  MHz,  and  ff  ~  60°,  we  calculate  the 
electron  density  n*  required  in  the  upper  hot  layer  if  the  emission  process  is  due  to  2nd  harmonic 
gyroemission  (n^  =  3.0  x  10“*  cm~^),  3rd  harmonic  gyroemission  (ne  =  6.8  x  10®  cm"®),  and  4th 
harmonic  gyroemission  (n^  =  1.2  x  10®  cm"®). 

A  second  possible  model  consists  of  two  layers  of  plasma:  the  lower  hot  layer  seen  by  XRP  and 
a  cooler  overlying  layer  where  gyroresonance  absorption  takes  place.  In  this  cool  absorbing  layer, 
neither  mode  is  optically  thick  and  since  the  X-mode  is  more  attenuated  than  the  0-mode,  the 
resulting  emission  is  dom.‘nated  by  0-mode.  The  observed  emission,  however,  seems  to  be  X-mode 
(it  is  left-circularly  polariz.’d,  with  an  inward-directed  underlying  photospheric  field).  Nevertheless, 
the  microwave  emission  can  undergo  a  polarization  inversion  upon  traversing  a  “quasi-transverse” 
(QT)  layer  in  the  coronal  magnetic  field  (see  Fig.  8b).  This  QT  layer  occurs  where  the  magnetic 
field  becomes  transverse  to  the  Dne  of  sight,  i.e.,  the  longitudinal  field  component  Bz  becomes  zero. 
Previous  studies  have  provided  evidence  for  such  polarization  inversion  (Kundu  et  al.  1977,  Webb 
et  al.  1983,  Alissandrakis  and  Kundu  1984).  The  observed  polarization  is  inverted  when  coupling 
between  the  polarization  modes  is  weak. 

Zheleznyakov  (1970),  using  the  WKB  method,  obtains  the  following  expression  for  the  coupling 
parameter,  given  here  in  the  notation  of  Bandiera  (1982): 

2  In  2  (mec)“*  a>“* 

gS  1  ^5  I’  VO 

where  mg  and  e  are  the  rest  mass  and  charge  on  an  electron,  c  is  the  speed  of  light,  w  is  the 
observation  angular  frequency,  Ng  is  the  electron  number  density,  B  is  the  magnetic  field  strength, 
and  d0/ds  is  the  gradient  of  the  angle  between  the  field  and  the  line  of  sight;  all  units  are  in  CGS. 
The  microwave  emission  undergoes  a  polarization  inversion  when  C  <  1,  or 

4.766  X  I  f  K  I.  (8) 

In  and  near  the  location  of  the  XRP  spectroscopic  scan,  the  Sakurai  code  yields  a  QT  layer  at 
a  height  of  ~75,000  km.  The  average  total  magnetic  field  strength  obtained  from  the  code  is  17 
G,  and  the  angular  gradient  along  the  line  of  sight  is  1.1  X  10"*°  rad/cm.  Inserting  these  values 
into  the  expression  for  the  coupling  parameter  and  using  the  highest  observing  frequency  in  the  20 
cm  waveband  (1.665  GHz),  we  find  that  the  polarization  inversion  will  take  place  if  the  electron 
density  in  the  QT  layer  exceeds  8.4  x  10®  cm"®.  The  model  plage  electron  density  of  Hildebrandt 
et  al.  (1987)  at  a  height  of  75,000  km  is  1.4  x  10®  cm"®,  and  the  model  electron  densities  in  all 
of  the  coronal  structures  in  Hildebrandt  et  al.  (1987)  and  Pagagiannis  and  Kogut  (1975)  exceed 
5.4  X  10"  cm"®  at  the  same  height.  Thus  it  is  reasonable  to  expect  polarization  inversion  to  occur 
for  the  20  cm  microwave  emission  from  the  plage:  the  20  cm  emission  e.\iting  the  cool  absorbing 
plasma,  dominated  by  the  0-mode  (RCP),  is  inverted  in  the  QT  layer  so  that  the  0-mode  emission 
becomes  LCP,  tis  observed. 

Plcisma  and  magnetic  field  parameters  in  the  cool  absorbing  layer  can  be  calculated.  The  optical 
depths  of  the  cool  absorbing  layer  can  be  estimated  by  solving 

T^.O  ^  exp(-r,''''0)  ^  exp(-r,'^'-0)]  (9) 
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for  where  are  the  observed  X-  and  0-  mode  brightness  temperatures  at  a  particular 
frequency  (prior  to  polarization  inversion),  is  the  brightness  temperature  due  to  thermal 

bremsstrahlung  predicted  from  the  X-ray  observations,  and  Tc  is  the  electron  temperature  of  the 
cool  absorbing  layer.  Solving  for  we  obtzdn 

=  ln[{T,  -  T|«^^)/(2’c  -  (10) 

For  the  average  observed  brightness  temperatures  in  the  20  cm  waveband,  we  obtain  =  1.44 
and  =  0.18  for  Tc  =  3  x  10®  K,  which  leads  to  0  ~  62°.  This  angle  is  comparable  to  the  angle 
obtained  from  the  extrapolated  field.  Using  r/  ~  1.44,  Tc  =  3  x  10®  K,  i/  ~  1500  MHz,  and 
6  ~  60°,  we  calculate  the  electron  density  ne  required  in  the  cool  absorbing  layer  if  the  absorption 
process  is  due  to  2nd  harmonic  gyroabsorption  (nj  =  5.5  x  10®  cm"^),  3rd  harmonic  gyroabsorption 
(ne  =  1.2  X 10®  cm~^),  and  4th  harmonic  gyroabsorption  (ne  =  2.1  x  10^^  cm"®).  The  2nd-harmonic 
density  is  considerably  less  than  canonical  low  coronal  density  values,  but  cannot  be  ruled  out.  The 
3rd-harmonic  density  is  consistent  with  canonical  low  coronal  density  values,  and  is  also  consistent 
with  the  density  requirement  for  polarization  inversion  in  the  QT  layer  (cdthough  the  QT  layer 
is  ~  70,000  km  higher  in  the  atmosphere).  The  4th-harmonic  density  is  so  high  that  plasma 
suppression  would  prohibit  the  propagation  of  20  cm  electromagnetic  radiation,  so  this  can  be 
ruled  out.  The  absorption  which  produces  the  observed  polarization  at  the  location  of  the  XRP 
spectroscopic  scan  most  likely  occurs  in  the  3rd  harmonic  gyroresonance  absorption  region. 

The  situation  is  somewhat  different  at  the  location  of  the  20  cm  intensity  peak  in  the  plage. 
Here,  ~  30,000  km  west  of  the  location  of  the  XRP  spectroscopic  scan,  the  predicted  and  the 
observed  20  cm  intensity  peaks  are  cospatial.  The  predicted  brightness  temperature  is  2.2  x  10®  K, 
while  the  observed  LCP  brightness  temperature  is  1.9  x  10®  K  and  the  observed  RCP  brightness 
temperature  is  0.9  x  10®  K.  The  maximum  predicted  6  cm  brightness  temperature  is  7x  10®  K,  while 
the  observed  is  <  4.2  X 10®  K.  This,  again,  requires  the  presence  of  cooler  absorbing  plasma  overlying 
the  hot  plasma  observed  by  XRP.  Here,  however,  the  coronal  potential  field  is  not  sufficient  for  the 
gyroemission  mechanism  to  operate.  The  photospheric  field  is  less  than  100  G,  and  between  heights 
of  5000  and  15,000  km  the  total  e.xtrapolated  field  ranges  from  ~  60  to  ~  90  G.-  This  suggests  quite 
strongly  that  the  extrapolated  field  is  substantially  less  than  the  actual  magnetic  field  in  the  corona. 
This  is  possibly  a  manifestation  of  instrumental  insensitivity  to  photospheric  longitudinaJ  magnetic 
fields  less  than  ~  100  G.  It  is  quite  possible  that  compact,  highly  intense  photospheric  magnetic 
field  regions  were  not  detected  due  to  inadequate  instrumental  resolution.  It  also  cannot  be  ruled 
out  that  electric  currents,  and  hence  non-potential  magnetic  fields,  are  present  in  the  corona. 

These  models  represent  the  two  simplest  ways  to  describe  and  interpret  the  simultaneous,  multi¬ 
waveband  data  taken  of  AR  4906.  Simpler  models  can,  of  course,  explain  the  microwave  images, 
the  X-ray  spectra,  or  the  magnetic  extrapolations  individually  or  even  explain  any  two  of  these, 
but  the  two  models  explained  in  detail  above  and  depicted  in  Figure  8  can  describe  the  entire  data 
set. 


4  Conclusions 


A  number  of  important  conclusions  concerning  the  plasma  and  magnetic  field  structure  of  an  active 
region,  as  well  as  the  observational  methods,  result  from  this  study: 
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(1)  The  temperature  in  the  corona  over  an  active  region  can  diminish  with  height;  this  occurs 
over  both  the  sunspot  and  the  plage  in  AR  4906  where  a  different  line  of  reasoning  leads  to  this 
conclusion  for  each  structure.  In  the  case  of  the  sunspot,  the  cospatial  microwave  emission  in  the  6 
and  20  cm  wavebands  have  substantially  different  brightness  temperatures.  The  emission  in  both 
wavebands  is  attributed  to  optically  thick  thermal  gyroemission  in  the  high  magnetic  field  of  the 
sunspot.  The  6  cm  emission  requires  higher  fields  than  the  20  cm  emission  and,  hence,  originates 
lower  in  the  sunspot  atmosphere.  Where  the  6  and  20  cm  emission  overlap,  the  6  cm  brightness 
temperature  is  2.5  x  10®  K  while  the  20  cm  brightness  temperature  is  1.3  x  10®  K.  Thus  the  sunspot 
atmospheric  temperature  diminishes  from  2.5  x  10®  K  at  the  height  of  the  6  cm  emission  to  1.3  x  10® 
K  at  the  height  of  the  20  cm  emission.  In  the  case  of  the  plage,  the  conclusion  is  beised  upon  the 
discrepancy  between  the  observed  RCP  emission  at  20  cm  and  the  brightness  temperature  predicted 
from  the  XRP  line  flux  ratios.  Specifically,  this  requires  the  presence  of  overlying  cool,  absorbing 
plasma  in  order  to  reduce  its  predicted  brightness  temperature  of  T^^^^  ~  2x  10®  K  to  its  observed 
value  of  Tg^  ~  0.7  x  10®  K.  The  overlying,  cool  plasma  must  have  a  temperature  <  5  x  10®  K. 

(2)  A  loop  or  loop-like  structure  connects  the  sunspot  with  a  nearby  region  of  opposite  polarity 
in  the  east.  The  peaks  of  the  6  and  20  cm  emission  lie  along  this  loop,  with  higher  frequency 
emission  originating  farther  down  the  leg  of  the  loop  on  the  side  of  the  sunspot.  Based  upon  the 
potential  field  extrapolation,  the  20  cm  emission  is  due  to  a  combination  of  2nd  and  3rd  harmonic 
gyroemission,  putting  the  magnetic  field  strength  in  this  region  in  the  range  160  -  300  G. 

(3)  Based  upon  a  comparison  of  the  6  cm  and  X-ray  observations  with  calculations  of  the 
emission  from  a  model  loop,  we  deduced  the  plasma  and  magnetic  field  parameters  at  the  location 
of  the  6  cm  source.  An  electron  temperature  T*  =  2.5  x  10®  K,  an  electron  density  Ne  ~  1  x  10® 
cm~^,  and  a  magnetic  field  strength  of  583  G  (3rd  harmonic  gyroemission  at  4900  .\IHz)  result. 
The  potential  field  extrapolation  yields  this  magnetic  field  strength  at  a  height  ~  5000  km  above 
the  photosphere. 

(4)  These  results  have  allowed  us  to  obtain  good  estimates  of  the  Alfven  speed  at  two  locations 
within  the  loop.  .\t  the  6  cm  source  near  the  sunspot,  an  Alfven  speed  ~  40.000  km/sec  is  obtained. 
Higher  in  the  loop,  at  the  location  of  the  20  cm  emission,  the  .A.lfven  speed  is  found  to  range  from 
10,000  km/sec  -  20,000  km/scc.  These  values  are  substantially  higher  than  the  1000  -  2000  km/sec 
typically  quoted  for  the  Alfven  speed  in  the  corona.  At  both  locations  the  electron  gyrofrequency 
exceeds  the  plasma  frequency. 

(5)  The  extrapolated  potential  sunspot  field  can  be  approximated  by  that  of  a  vertically  oriented 
point  dipole  for  locations  in  which  B  ~  300  G.  For  locations  in  which  B  <  300  G,  nearby  pockets 
of  magnetic  field  perturb  the  sunspot  field  from  that  of  a  point  dipole. 

(6)  The  electron  temperature  and/or  density  are  lower  over  the  western  portion  of  the  sunspot 
(away  from  the  center  of  the  active  region)  than  they  arc  over  the  eastern  portion  of  the  sunspot. 

(7)  In  the  plage,  the  ratios  of  X-ray  emission  lines  are  used  to  calculate  the  electron  temperature 
and  column  emission  meaisure.  The  most  accurate  determination  of  these  parameters  is  at  the 
location  of  the  Fe  XVII  peak  flux,  where  spectroscopic  scan  data  were  obtained.  The  Fe  XVII 
emissivity  needs  to  be  divided  by  a  correction  factor  ~  2.75.  This  is  somewhat  higher  than  but 
comparable  to  the  ~  2.25  factor  obtained  in  CoMSlOC  I. 
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(8)  The  predicted  20  cm  brightness  temperature  in  the  plage  is  reasonably  consistent  with 
the  observed  for  the  LCP  emission  but  is  much  higher  for  the  RCP  emission.  We  have  proposed 
two  possible  models  to  explain  this  discrepancy.  Both  models  are  as  simple  as  possible  yet  still 
incorporate  all  three  sets  of  data  for  AR  4906:  The  first  model  has  a  low  layer  of  hot  plasma 
(where  the  X-ray  emission  originates)  covered  by  an  optically  thick  cool,  absorbing  layer.  On  top 
of  this  cool  plasma  is  a  second  hot  layer  which  is  optically  thin  to  thermal  bremsstrahlung  but 
emits  strong,  polarized  microwave  emission  since  it  is  embedded  in  a  magnetic  field  sufficient  for 
gyroresonance  emission  at  20  cm.  Since  the  longitudinal  component  of  the  ma^etic  field  inverts 
at  a  height  of  ~75,000  km,  the  density  at  this  height  cannot  exceed  ~  8  x  10®  cm“®;  otherwise, 
the  polarization  will  invert  in  this  QT  layer  and  the  sense  of  polarization  will  be  opposite  that 
observed.  (It  should  be  pointed  out  that  the  model  coronal  densities  of  Hildebrandt  et  al.,  1987, 
and  Papagiannis  and  Kogut,  1975,  greatly  exceed  8  x  10®  cm“^  at  a  height  of  75,000  km.)  The 
second  model  has  a  layer  of  cool,  but  not  fully  thick,  absorbing  plasma  overlying  the  hot  plasma 
from  which  the  predicted  microwave  emission  originates.  This  cool  plasma,  embedded  in  a  magnetic 
field  sufficient  for  gyroresonance  absorption,  attenuates  the  X-mode  more  strongly  than  the  0- 
mode.  The  polarization  is  subsequently  inverted  in  the  overlying  QT  layer  in  the  magnetic  field. 
Appropriate  coronal  magnetic  field  strengths,  as  well  as  the  QT  layer,  are  extrapolated  from  the 
photospheric  longitudinal  magnetogram  with  the  Sakurai  (1982)  code.  A  similar  result  has  been 
obtained  from  the  analysis  of  the  CoMStOC  II  data  (Schmelz  et  al.  1991). 

(9)  Elsewhere  in  the  plage,  at  the  observed  20  cm  intensity  peak,  we  find  that  the  magnetic  field 
required  to  explain  the  observations  is  higher  than  the  coronal  potential  field  extrapolated  from 
the  photospheric  longitudinal  magnetogram.  A  similar  discrepancy  has  been  found  in  the  analysis 
of  other  CoMStOC  data. 
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TABLE  1 

Plasma  Parameters  from  XRP  Spectroscopic  Scan  Data 


Line  Ratio 

log(EM)  (cm“®)(^^ 

0  VIII/Mg  XI 

3.01  ±  0.16 

28.85  ±0.08 

0  Vin/Fe  XVIII 

2.83  ±  0.19 

28.87  ±0.09 

Ne  IX/Fe  XVHI 

2.94  ±  0.22 

28.85  ±0.09 

Mg  XI/Fe  XVm 

2.61  ±  0.50 

28.89  ±0.24 

Fe  XVII/Fe  XVHI 

3.12  ±  0.27 

28.83  ±0.06 

weighted  mean^^) 

2.95  ±  0.20 

28.85  ±0.19 

The  quoted  uncertainties  on  the  individual  values  of  Te  reflect  counting  statistics  only. 

The  column  emission  measure  is  computed  from  the  0  VIII  flux  only;  the  quoted  uncertainties 
on  the  individual  values  reflect  counting  statistics  only. 

The  quoted  uncertainties  on  the  weighted  means  of  and  EM  are  the  root  mean  squares 
of  the  statistical  and  systematic  uncertainties.  For  temperature,  the  systematic  uncertainty  is  esti¬ 
mated  as  the  standard  deviation  of  the  individual  values.  For  the  emission  measure,  the  systematic 
uncertainty  is  estimated  as  the  standard  deviation  of  the  emission  measures  calculated  for  the  0 
VIII,  Ne  IX,  Mg  XI,  and  Fe  XVIII  fluxes  at  the  mean  Tg. 
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TABLE  2 

Predicted  Microwave  Intensities  from  XRP  Spectroscopic  Scan  Data  And  Observed  Microwave 

Intensities  from  the  VLA 


Freq.  (MHz) 

IPRED  (MK) 

lOBS 

Tgl^  (MK) 

Tgf/  (MK) 

1375 

2.17  ±0.49 

1.2 

1.5 

0.9 

1420 

2.10  ±0.49 

1.1 

1.5 

0.7 

1550 

1.91  ±  0.49 

1.4 

2.0 

0.8 

1665 

1.75  ±  0.49 

1.1 

1.7 

0.5 

4600 

0.31 

<0.42 

— 

— 

4700 

0.30 

<0.42 

— 

— 
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Figure  Captions 

Figure  1.  Observations  of  the  CoMStOC  IV  active  region,  AR  4906.  Celestial  north  is  up, 
and  west  is  to  the  right.  All  eight  boxes  cover  the  same  field  of  view,  for  ease  in  coalignment.  The 
tick  speuung  is  one  minute  of  arc,  the  crosses  mark  the  locations  of  the  phase  center  of  the  VLA 
observations,  and  “X”  marks  the  location  of  the  XRP  spectroscopic  scan,  (a)  Frequency-averaged 
map  of  I  (intensity)  at  20  cm.  The  contour  levels  are  0.21,  0.41,  0.62,  0.83,  1.03, 1.24, 1.45,  1.66, 
1.86,  cind  2.07  xlO®  K.  (b)  Frequency-averaged  map  of  I  at  6  cm.  The  contour  levels  are  0.42, 
0.63,  0.84,  1.05,  1.26,  1.47,  1.68,  1.89,  and  2.10  xlO®  K.  (c)  XRP  soft  X-ray  map  of  AR  4906  in 
the  bright  Fe  XVII  line  at  15.01  A.  The  image  was  accumulated  during  a  A'  square  raster  with 
10"  steps  which  took  10  min  to  complete.  The  map  has  been  smoothed  with  a  3  x  3  boxcar  to 
improve  statistics.  Contour  levels  are  5,  10,  15,  20,  25,  30,  35,  and  40  counts/s.  (d)  Photospheric 
lon^tudinal  m^netogram  obtained  from  the  Beijing  Observatory,  with  superimposed  magnetic 
field  lines  calculated  with  the  potential  field  extrapolation  code  of  Sakurai.  The  contour  levels  are 
200  and  500  G,  where  the  dashed  contours  represent  inward  directed  fields  and  the  solid  contours 
represent  outward  directed  fields,  (e)  Electron  temperature  and  (f)  log  of  the  column  emission 
measure  calculated  from  line  flux  ratios  of  emission  lines  observed  with  XRP.  The  contours  in  (e) 
are  2.2, 2.4, 2.6,  2.8,  and  3.0  x  10®  K;  the  contours  in  (f )  are  28.7, 28.9,  and  29.1.  The  temperature 
and  emission  measure  di<^ostics  were  available  only  in  the  plage,  (g)  Brightness  temperature 
map  at  20  cm  calculated  for  the  thermal  bremsstradilung  emission  mechanism  using  the  electron 
temperature  and  column  emission  measure  maps  in  (e)  and  (f).  The  contours  are  1.7,  1.9,  and 
2.1  xlO®  K.  (h)  Brightness  temperature  map  at  6  cm  calculated  for  the  thermal  bremsstrahlung 
emission  mechanism  using  the  electron  temperature  and  column  emission  measure  maps  in  (e)  and 
(f).  The  contours  are  0.1,  0.3,  and  0.5  xlO®  K. 

Figure  2.  I  and  V  maps  at  each  of  the  four  frequencies  in  the  20  cm  waveband.  The  I  contour 
levels  are  0.53,  0.79, 1.05, 1.31,  1.57, 1.83, 2.09,  and  2.35  xlO®  K  for  all  frequencies.  The  V  contour 
levels  are  0.29  and  0.38  xlO®  K  at  1375  MHz;  0.27,  0.36,  and  0.46  xlO®  K  at  1420  MHz;  0.31, 
0.41,  0.51,  and  0.61  xlO®  K  at  1550  MHz;  and  0.40,  0.53,  and  0.66  xlO®  K  at  1665  MHz.  Celestial 
north  is  up,  and  west  is  to  the  right.  The  tick  spacing  is  one  minute  of  arc,  and  the  crosses  mark 
the  location  of  the  phase  center  of  the  VLA  observations. 

Figure  3.  Sequence  of  contour  maps  showing  the  extrapolated  magnetic  field  strength  at 
heights  of  (a)  0  km,  (b)  5000  km,  (c)  10000  km,  and  (d)  15000  km  above  the  photosphere. 
Celestial  north  is  up,  and  west  is  to  the  right.  The  numbers  along  the  axes  are  distances  from  solar 
disk  center  in  units  of  10^  km  in  the  plane  of  the  sky.  The  contour  levels  in  (a)  and  (b)  are  150, 
250,  350,  500,  and  800  G.  The  contour  levels  in  (c)  and  (d)  are  150,  200,  250,  350,  and  500  G. 
Crosshmrs  mark  the  phase  centers  of  the  VLA  6  and  20  cm  maps,  to  assist  in  coaligning  images. 
A  40"  line  segment  is  shown  in  (d)  to  assist  in  scaling  and  coaligning. 

Figure  4.  Magnetic  field  of  a  model  sunspot  consisting  of  a  vertically  oriented  point  dipole 
buried  below  the  photosphere.  Numbers  along  the  x  and  y  axes  represent  the  distance  from  the 
point  dipole  in  units  of  10^  km.  The  x-axis  is  in  the  plane  of  the  photosphere,  and  the  y-ads 
is  perpendicular  to  that  plane.  Solid  lines  are  field  lines,  dashed  lines  numbered  1  through  5  are 
contours  of  constant  magnetic  field  strength,  and  the  dotted  horizontal  lines  represent  heights  above 
the  photosphere  of  5000, 10000,  and  15000  km.  In  (a)  the  dashed  lines  are  ma^etic  field  contours 
corresponding  to  the  1st  through  5th  harmonics  of  1665  MHz:  594,  297,  198,  149,  and  119  G.  In 
(b)  the  dashed  lines  are  magnetic  field  contours  corresponding  to  the  1st  through  5th  harmonics 
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of  4900  MHz:  1749, 875, 583, 437,  and  350  G. 


Figure  5.  I  and  V  maps  calculated  for  1665  MHz  and  4900  MHz  for  the  sunspot  model  whose 
mj^etic  field  is  shown  in  Figure  4.  Numbers  along  the  x  and  y  axes  are  in  units  of  arcsec,  and  the 
crosshmrs  intersect  at  the  center  of  the  sunspot.  Contour  line  thickness  increcises  with  increasing 
contour  value.  In  (a)  the  1665  MHz  I  contours  are  0.84,  1.37,  and  1.90  xlO®  K,  and  in  (b)  the 
1665  MHz  V  contours  are  0.27,  0.53,  and  0.80  xlO®  K.  In  (c)  the  4900  MHz  I  contours  are  0.66, 
1.28,  and  1.90  xlO®  K,  and  in  (d)  the  4900  MHz  V  contours  are  0.31,  0.62,  and  0.93  xlO®  K. 

Figure  6.  Portion  of  the  magnetogram  showing  the  locations  of  the  6  and  20  cm  emission  along 
the  field  lines  which  outline  the  loop-like  structure  connecting  the  sunspot  with  the  nearby  region  of 
opposite  polarity.  Celestial  north  is  toward  the  top,  and  west  is  toward  the  right.  .4  segment  of  10" 
length  is  shown  for  scaling  purposes.  Shaded  areas  show  where  the  photospheric  longitudinal  field 
exceeds  200  G.  The  solid  magnetic  field  contours  represent  outward  directed  fields  of  200  and  500 
G,  while  the  short-dashed  ms^etic  field  contours  represent  inward  directed  fields  of  200  «ind  500 
G.  The  sunspot  is  the  large  shaded  area  in  the  upper  right.  The  long-dashed  contour  corresponds 
to  an  intensity  I  of  0.79  x  10®  K  at  the  representative  20  cm  microwave  frequency  of  1375  MHz. 
The  “-f”  symbols  mark  the  locations  of  the  peak  brightness  temperature  of  each  of  the  four  closely- 
spaced  frequencies  in  the  20  cm  waveband.  The  location  of  the  peak  brightness  temperature  shifts 
systematicedly  closer  to  the  footpoint  inside  the  sunspot  with  increasing  frequency:  the  “-f”  fartliest 
to  the  left  corresponds  to  the  1375  MHz  emission,  and  the  "-f-"’  farthest  to  the  right  corresponds 
to  the  1665  MHz  emission.  The  dotted  contour  corresponds  to  an  intensity  I  of  1.56  x  10®  K  at 
4900  MHz  in  the  6  cm  waveband.  The  dot-dash  contours  correspond  to  a  V  of  0.39  x  10®  K  at  4900 
MHz. 

Figure  7.  Theoretical  4900  MHz  I  and  V  maps  for  the  model  loop  described  in  the  text.  The 
height  of  the  loop  apex  (16000  km),  the  footpoint  separation  (50000  km),  and  the  loop  thickness 
at  the  apex  (7700  km)  were  all  estimated  from  the  extrapolated  potential  field  lines.  The  cross 
section  of  the  loop  is  taken  to  be  circular  at  the  apex.  The  model  loop  has  a  constant  temperature 
of  2.5  X  10®  K,  a  density  of  1  x  10^  cm"^,  and  a  field  strength  of  583  G  (3rd  harmonic)  in  the 
footpoints  at  a  height  of  ~  -5000  km  above  the  photosphere.  The  I  contours  are  0.52.  l.O  l.  1.56. 
and  2.34  x  10®  K.  The  V  contours  are  0.39. 0.52, 0.64,  and  1.00  x  10®  K.  The  model  I  and  V  sources 
on  the  right  side  of  the  figure  reproduce  the  observed  6  cm  I  and  V  sources  fairly  well. 

Figure  8.  Illustration  of  the  two  models  used  to  explain  the  observed  polarization  of  the 
20  cm  microwave  emission  from  the  plage.  The  arrows  labeled  with  an  “R”  correspond  to  RCP 
emission  and  the  arrows  labeled  with  an  “L”  correspond  to  LCP  emission.  The  propagation  mode 
is  indicated  by  an  “X”  or  an  “0”  in  parentheses.  The  propagation  mode  is  determined  by  the 
orientation  of  the  local  ms^etic  field,  indicated  by  arrows  labeled  B-.  The  lengths  of  the  arrows 
labeled  R  and  L  indicate  the  relative  brightness  temperatures  of  the  emergent  emission  in  the 
given  mode  of  polarization,  (a)  The  hot/cool/hot  layer  model.  Emission  with  high  brightness 
temperature  and  low  polarization  (Tg ~  2  x  10®  K)  emerges  from  the  lowest  layer.  This  emission 
is  then  absorbed  by  the  overlying  cooler  layer,  from  which  low  brightness  temperature  emission 
{Tg'^  <  5  X  10®  K)  emerges.  Highly  polarized  emission  {T^  ~  2  x  10®  K,  ~  7  x  10®  K) 
emerges  from  the  highest  layer,  which  is  a  hot  plasma  embedded  in  a  magnetic  field  which  is 
sufficiently  strong  for  gj'roresonance  emission,  (b)  The  hot/cool/QT  layer  model.  Emission  with 
high  brightness  temperature  and  low  polarization  {Tg'^  ~  2  x  10®  K)  emerges  from  the  lowest  layer. 
This  emission  is  then  attenuated  in  the  overlying  cooler  layer,  in  which  magnetic  fields  sufficient 
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for  gyroresonance  absorption  cause  the  optical  depth  (and  hence  the  attenuation)  to  be  greater  for 
the  X-mode  (L)  than  for  the  0-mode  (R).  This  results  in  the  emission  of  0-mode-dominated  (R) 
emission.  The  polarization  is  then  inverted  to  its  observed  (L)  value  upon  traversal  of  the  overlying 
quasi-transverse  (QT)  layer  in  the  magnetic  field  overlying  the  plage. 
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